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Abstract 

Correlations between charged particles in deep inelastic e^p scattering have been studied in 
the Breit frame with the ZEUS detector at HERA using an integrated luminosity of 6.4 pb""'^. 
Short-range correlations are analysed in terms of the angular separation between current-region 
particles within a cone centred around the virtual photon axis. Long-range correlations between 
the current and target regions have also been measured. The data support predictions for the 
scaling behaviour of the angular correlations at high and for anti-correlations between the 
current and target regions over a large range in and in the Bjorken scaling variable x. Analytic 
QCD calculations and Monte Carlo models correctly describe the trends of the data at high Q^, 
but show quantitative discrepancies. The data show differences between the correlations in deep 
inelastic scattering and e"'"e~ annihilation. 
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1 Introduction 



This paper reports a study of both short-range and long-range correlations in neutral-current deep 
inelastic e^p scattering events measured at a centre-of-mass energy of -y/s = 300 GeV with the 
ZEUS detector at HERA. 

In deep inelastic scattering (DIS) at HERA energies, final-state charged particles have been stud- 
ied in terms of fragmentation functions transverse momentum spectra |6| and particle 
correlations 0-^]- The short-range correlations between final-state particles at small phase-space 



separations are sensitive to the structure of many-particle inclusive densities [10|. The study of 
long-range correlations can aid the understanding of the interdependence between distant phase- 
space regions; these measurements complement and extend the studies of the global multiplicity 
distributions that have already been reported in DIS QpllTif at HERA energies. 



Two-particle angular correlations have recently been studied in e^e collisions by the DELPHI 
Collaboration [12| at LEP and many-particle correlations have also been measured by the L3 [^] 



and OPAL Collaborations |14|. Such studies provide more detailed insight into the QCD models 



(see reviews [15,1^) based on the hypothesis of Local Parton-Hadron Duality, which states that the 



single-particle distributions of hadrons are proportional to the corresponding parton spectra |17] 



The analytic predictions, derived in the Double Log Approximation (DLA) |jT8-2C], have shown 
discrepancies with the observed angular correlations as well as with the fluctuation measurements 
]l2|,p!3[. The DLA calculations take into account the contributions from the infrared and collinear 
singularities of the gluon emissions and the angular ordering (colour coherence), but neglect energy- 
momentum conservation in gluon splittings and qq production as well as flnite-energy effects in the 
QCD parton cascade. However, Monte Carlo (MC) models which take these effects into account 
and which, in addition, contain hadronisation followed by resonance decays, also show discrepancies 



with the data |13, 1 



DIS provides a unique possibility to confront both the analytic results and MC models with data 
at various energy scales. The main motivation for such studies is to investigate the applicability 
of the perturbative analytic calculations to DIS data at HERA energies. The two-particle angular 
correlations in the current region of the Breit frame are also compared with the correlations in 
a single hemisphere of e^e~ annihilation in order to study a possible universality of two-particle 
spectra in these processes. Recent comparisons of single-particle spectra have shown considerable 
similarity between DIS at high and e^e~ annihilation [Q-Q,^. 

Correlations between distant phase-space regions provide important information on multi-hadron 
production that cannot be studied by measuring observables separately in each phase-space region. 
In DIS, a measurement of the correlations between the hemisphere populated by hadrons originating 
from the struck quark and that containing predominantly particles from the proton remnant can 
be made in a way analogous to "forward-backward" correlations. Such long-range correlations in 
rapidity have been studied for many years [|, 21~p^. Results for e'^e~ and /i+p processes indicated 



that such correlations are small at low energies |21|, |22(|. At LEPl energies, positive long-range 
correlations were observed, mainly due to heavy quark pair production p^-p5[. For vp and vp 



processes, the correlations are small and negative For pp, |27| and vr p, K~^p |£8[ collisions, 
the correlations are positive and increase with ^/s. Recently, the HI Collaboration has shown that 
the forward- backward correlations measured in the 7*-pomeron centre-of-mass system of diffractive 
e'^p collisions are positive Q. In contrast, in the Breit frame in DIS, negative long-range correlations 
due to the kinematics of first-order QCD effects have been predicted [pSj]. 



1 



2 Definitions and analytic QCD results 



2.1 DIS kinematics in the Breit frame 

The event kinematics of DIS processes are determined by the negative squared 4-momentum transfer 
at the lepton vertex, = —q^ = — (/c — /c )^ {k and k denote the 4-momenta of the initial and 
final-state positron, respectively), and the Bjorken scaling variable x = Q'^/{2Pp ■ q), where Pp is 
the 4-momentum of the proton. The fractional energy transfer to the proton in its rest frame, y, 
is related to these two variables by y ~ Q'^/xs, where -^/i is the positron-proton centre-of-mass 
energy. 

The correlations are studied in the Breit frame which provides a natural system to separate the 
radiation from the outgoing struck quark from that associated with the rest of the hadronic final 
state. In this frame, the exchanged virtual boson is completely space-like and has a momentum 
Q = {Qo,<lx,<lY,qz) = (0, 0,0, — Q) as shown in Fig. ||. For the Quark-Parton Model (QPM), the 
incident quark carries momentum p^^*^^^ = Q/2 in the positive Z-direction and the outgoing 
struck quark carries p^^'^^^ = —Q/2 in the negative Z-direction. The phase space of the event 
can be divided into two regions. All particles with negative p^^^^^ form the current region. These 
particles are produced by the fragmentation of the struck quark, so that the current region is 
analogous to a single hemisphere of e^e~ annihilation. Particles with positive p^''^'* are assigned 
to the target region associated with the incoming proton. The longitudinal phase space in p^^*^^^ 
available for the target-region particles extends to p^'^^^^ = (1 — x)Q/2x. This means that the 
current and target regions are highly asymmetric at small x. Thus the Breit frame differs from the 
hadronic centre-of-mass (7*p) frame and the centre-of-mass frame of e~^e~ annihilation. 

This simple picture of DIS in the Breit frame is changed when leading-order QCD processes, Boson- 
Gluon Fusion (BGF) and QCD-Compton (QCDC) scattering, are considered. In such processes, 
the collision in the Breit frame is no longer collinear, although the current and target regions are 
still well-defined operationally pl| . 



2.2 Two-particle correlations 

The correlations between charged particles in the current region are studied using the angle between 
any two particles (see Fig. |l]) . For each pair of particles inside a cone of half-angle G centred around 
the Breit frame axis, the relative angle 9i2 is determined. For a given QCD scale, A, this angle is 
transformed into the variable 

ln(P sine/A)' 

For DIS in the Breit frame, P = Q/2 is chosen, corresponding to the outgoing quark momentum 
in the QPM. Note that A is an effective scale parameter which is not related to Aj^ [|T^. The 
two-particle inclusive density 

is determined, where N^v is the number of events and ((inpair/de)fe is the number of particle pairs 
in the interval e to e + 6e. The definition (||) is used since the DLA predictions can be naturally 
expressed in this scaling variable |18,2^ and the logarithmic dependence of (|l|) expands small values 



of 9i2. Note that our definition of e is different from the scaling variable with the denominator 



2 



ln(P0/A) used previously [|l2,13|: for large 0, the definition (Q) is better suited to compare the 



DLA predictions with the data |32|. 

In this paper we measure the relative angular distribution f(e) and the correlation function r{e): 

The variable f(e) is p(e) normalised using the average charged multiplicity (n(B)) in the cone. Even 
for uncorrelated particle production, f(e) depends on e in a manner which is determined by the 
single-particle spectra. This dependence is reduced in r(e), which is normalised instead by Pmixi^), 
calculated using particles from diffcTCTit events. Thus Pmix 

(e) does not contain any of the dynamical 
correlation present among particles from the same event, but is sensitive to effects determined by 
the single-particle spectra. To obtain Pmmi^), each of the original tracks in an event is replaced by 
a track selected at random from all the other events. This is performed after the transformation of 
the real events to the Breit frame. These "fake" events are then used to calculate Pmixi^) in exactly 
the same way as p{e) is calculated from the real events. It should be noted that each fake event 
has, by construction, the same multiplicity (n) as the real event, but there is no requirement that 
the fake tracks are drawn from events of multiplicity close to n. This is the same procedure as used 



by DELPHI [12| and allows direct comparison with theoretical calculations. However, it has the 
consequence that the value of r(e) is influenced by the mixing of events with different multiplicities 
in Pmix(e)- For this reason many correlation studies, particularly in hadronic collisions, have used 
the "semi-inclusive" correlation [10|, r(")(e), for which both p{e) and Pmixi^) are evaluated at a 
fixed multiplicity n. Monte Carlo studies indicate that, for the data studied in this paper, the 
difference between r(e) and r^^\e) is quite pronounced for low multiplicity ((n(G)) < 2), but 
largely disappears at high multiplicity. 

Analytic calculations for f{e) and r{e), performed in the DLA at asymptotic energies, can be 
expressed in the form ||T8| : 

Y.--^^lil=. 2,3(1 -<>.Me)), (4) 
2Y^ln(Psm0/A) 

Inr(e) 



Y ^ ] ' ^ 2/3(u;(e) - 2^/1^), (5) 

y'[n[F sm fc)/A) 



where 



uj{e) = 2Vl -e (1 -ln(l -e)/8) , (6) 

and (3 and b are given by 

/? = 6(117Vc-2n/)-^/^ 6 = 2/?Vln(-Psine/A) (7) 
with Nc = 3. The effective number of flavours, nj, is chosen to be three since the main contribution 



to the current-region parton cascade comes from light quarks [ 33 1 . These values of Nc and n f give 
P = 1.15. 

The analytic QCD calculations predict a scaling behaviour, i.e. Y and Y are functions of e only 
and depend neither on the initial parton momentum P = Q/2 nor on the angle Q. According to 
the calculations, both Y{£) and Y(e) rise with increasing e. The increase of Y{e) determines the 
behaviour of f(e) which is a strongly decreasing function due to the decrease in the number of 
partons (and, hence, parton pairs) with increasing e (decreasing ^12). In contrast, r(e) is predicted 



3 



to rise with increasing e, reflecting an increased sensitivity of r(e) to two-particle correlations. At 
large particle separation, i.e. 6*12 > 6 (e < 0), there are no two-particle correlations, r{e < 0) = 1 



2.3 Current-Target correlations 

A simple way to measure the interdependence between the number of charged particles in the 
current and target regions, Uq and rit, is to use the correlation coefficient k: 

K = a~'^a~^cov{nc,nt), cov(nc, nt) = (ncTi-t) - (nc)(nt), (8) 

where cic and at are the standard deviations of the multiplicity distributions in the current and tar- 
get regions, respectively. For positive correlations, k, is positive; for anti-correlations it is negative. 
The advantages of using definition (|8|) lie in the simplicity of the boundary conditions, —1 < k < 1, 
which allows a quantitative estimate of the strength of the correlations, and its low sensitivity to 
the losses of particles from the proton remnants, which often escape undetected down the beam 
pipe. The latter property comes from the fact that cov(nc, rit) and the product cXc (Tt have a similar 
dependence on the total number of measured tracks, so that their ratio has only a small dependence 
on the average track multiplicity in the target region (see Sect. 



If hadronisation effects are neglected, the QCDC and BGF processes lead to [29| 

cov(ne,nt) ~ -AiiQ^)R,{Q^x)-A2{Q^)R2{Q^x), (9) 

where Ai{Q'^) and A2{Q'^) are x-independent functions determined by the multiplicity of the outgo- 
ing partons. Ri{Q'^,x) is the probability for back-to-back jet events (i.e. with one jet in the current 
and one in the target region) and R2{Q'^,x) is the probability for an event without jet activity in 
the current region. 



At small Q^, the values of Ai{Q'^) and ^2(Q^) are positive [^. This leads to a negative covariance 
(P) and thus to anti-correlations between the current and the target regions, which can be measured 
using eq. (^. At asymptotically high Q^, the value of Ai{Q'^) can be negative. In this case the 
current-target correlations are not strong and can even be positive. 

Note that at small fixed Q^, cov(nc,nt) is sensitive to the BGF rate as a function of x which, in 



turn, depends on the gluon density inside the proton |2£ 



3 Experimental setup 

The data were taken in 1995 at the positron-proton collider HERA with the ZEUS detector. During 
this period, the energy of the positron beam was = 27.5 GeV and that of the proton beam was 
820 GeV. The integrated luminosity used for the present study is 6.4 pb~^. 

ZEUS is a multi-purpose detector described in detail in |3^. Of particular importance in the 
present study are the central tracking detector and the calorimeter. 

^This statement follows from the analytic DLA calculation, which takes into account the angular ordering but 
neglects momentum conservation in the gluon splitting. On the other hand, Monte Carlo models which incorporate 
momentum conservation exhibit long-range two-particle correlations in this phase-space region ||l8|. 
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The central tracking detector (CTD) is a cylindrical drift chamber with nine superlayers covering 
the polar angleQ region 15° < 9 < 164° and the radial range 18.2-79.4 cm. Each superlayer consists 
of eight sense wire layers. The transverse momentum resolution for charged tracks traversing all 
CTD layers is (t{pi_)/p± = 0.0058p_L © 0.0065 © 0.0014/p_L with p±_ being the track transverse 
momentum (in GeV). The single hit efficiency of the CTD is greater than 95%. 

The CTD is surrounded by the uranium-scintillator calorimeter which is divided into three parts: 
forward (FCAL), barrel (BCAL) and rear (RCAL). The calorimeter is longitudinally segmented 
into electromagnetic (EMC) and hadronic (HAC) sections. The energy resolution of the calorimeter 
under test beam conditions is (Te/E = 0.18/\/^ for electrons and ge/E = 0.35/\/^ for hadrons 
(with E in GeV). 



4 Data selection 

The kinematic variables x and can be reconstructed using a variety of methods. In order to 
determine these variables, we use the measurement of energy and angle of the scattered positron, 
the double angle and the Jacquet-Blondel methods [^]. These three methods are further denoted 
by the subscripts e, DA and JB, respectively. In the double angle method, the variables x, and 
y are reconstructed using the angles of the scattered positron and the hadronic energy flow. In the 
Jacquet-Blondel method, used only in the event selection, the kinematic variables are determined 
entirely from the hadronic system. The boost vector from the laboratory to the Breit frame is 
determined using the double angle method, which is less sensitive to systematic uncertainties in 
the energy measurement than the other methods. In the reconstruction of the Breit frame all 
charged particles are assumed to have pion masses. 



The triggering and online event selections are identical to those used in [36|. To select neutral- 
current DIS events the following cuts are applied: 



• > 10 GeV, E^ being the energy of the scattered positron. 

• Qda > 10 GeV2. 

• 35<(5 = ^ii^i(l — cos 9i) < 60 GeV, where Ei is the energy of the i^^ calorimeter cell and 6i 
is its polar angle with respect to the beam axis. 

• ye< 0.95. 

• yjB > 0.04. 

• Z-component of the primary vertex position determined from the tracks fitted to the vertex 
is in the range —40 < .^vertex < 50 cm. 

• A timing cut requiring that the event time measured by the RCAL is consistent with an e~^p 
interaction. 

^ZEUS has a right-handed coordinate system in which X = Y = Z — is the nominal interaction point. The 
positive Z-axis is along the direction of the proton beam. The X-axis is horizontal, pointing towards the centre of 
HERA. The polar angle 9 is defined with respect to the positive ^-direction. 
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• The impact point {X, Y) of the scattered positron in the RCAL has to lie outside a square 
of 32 X 32 cm^, centred on the beam axis. 

About 350k events satisfy the above cuts. Tables || and ^ give the numbers of selected events for 
the kinematic regions used in this analysis. 

The present study is based on the CTD tracks fitted to the vertex. The scattered positron was 
removed from the track sample. In addition, the following cuts to select tracks were imposed: 

• Transverse momentum p± > 150 MeV. 

• I r/ |< 1.75, where r] is the pseudorapidity given by — ln(tan(0/2)) with 6 being the polar 
angle of the charged track. 

These cuts restrict our study to a well-understood, high-acceptance region of the CTD. 

To understand the uncertainties in our results, a subsample of DIS events containing a single 
jet was selected using calorimeter information. The cone algorithm |^] was used with radius 
R = drf' + Scf)^ = 1, where 5ri and 5(f) are the differences of pseudorapidity and azimuthal angles 
of energy deposits in the calorimeter with respect to the jet direction. Events with a single jet 
(in addition to particles in the proton remnant) are selected if the jet transverse energy is 
larger than 4 GeV. The jet pseudorapidity is restricted to the region | r/jet |< 2, where jets are well 
reconstructed. After the transformation to the Breit frame, only jets which belong to the current 
region of the Breit frame are considered. 



5 Event simulation 



To determine the corrections for selection and acceptance losses, event migration between 

bins, and track migration between the current and target regions due to mis-reconstruction of the 

Breit frame, a sample of neutral-current DIS events is generated with ARIADNE 4.8 using 



tuned parameters |^. This MC is based on the colour-dipole QCD model supplemented with the 
Boson-Gluon Fusion process. The hadronisation is described by the JETSET model [^^. Hadrons 
with lifetime cr > 1 cm are treated as stable. The GRV94 HO [41| parameterization of the 
proton parton distribution functions is used. The MC events obtained by this procedure define the 
generator-level sample. 

For the detector-level sample, the events are generated with ARIADNE using the DJANGO 6.24 [H^] 



program based on HERACLES 4.5.2 [43| in order to incorporate first-order electroweak corrections. 



The events are then processed through a simulation of the detector using GEANT 3.13 |44] to take 
into account particle interactions with the detector material, particle decays, event and track mi- 
grations, double counting of single tracks, resolution, acceptance of the detector and event selection. 
The detector-level MC events are processed with the same reconstruction program as used for the 
real data. 



In addition to ARIADNE, the LEPTO 6.5 g§ and HERWIG 5.9 |4|] generators with tuned pa- 



rameters [39 1 are also used to compare with the data. The parton cascade in LEPTO is based 
on the matrix element calculation matched to a parton shower (MEPS) according to the DGLAP 
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equation. LEPTO orders the parton emissions in invariant mass with an additional angular con- 
straint to ensure coherence. As in ARIADNE, the hadronisation in LEPTO is described by the 
JETSET model. HERWIG also has a parton shower based on the DGLAP equation, but parton 
emissions are ordered in angle. The hadronisation is described by a cluster model. 

One of the sources of two-particle correlations is the Bose-Einstein interference between identical 
particles. By default, this effect is turned off in the JETSET model. The simulation of the 



Bose-Einstein effect is absent in HERWIG. According to our studies and the DELPHI results |12], 
the Bose-Einstein interference has a small effect (less than 2% of relative change) on the angular 
correlations. 



6 Resolution and acceptance 

For the present analysis, it is important to understand the two-particle resolution of the CTD. The 
resolution is determined in the Breit frame as the variance of the absolute difference between 612 
measured at the generator-level and the detector-level Monte Carlo defined above. It is found that 
to ensure that the two particles are resolved, the angle 612 between two tracks in the current region 
of the Breit frame should not be smaller than two degrees. For a given O, P and A, this limit 
determines the maximum values of e used in this analysis. 

For the cuts used, the event and current-region track acceptances are both in the range 70 — 
90%, depending on the {Q'^,x) region studied. This is sufficiently good to measure the angular 
correlations reliably. 

For the target region, the track acceptance lies in the range 20 — 30%. The low acceptance re- 
sults mainly from the | |< 1.75 restriction. Using the ARIADNE model, the low target-region 
acceptance is found to have a small effect on the current-target correlations measured according to 
(^): the restriction | r] |< 1.75 decreases the absolute value of cov(nc,nt) by 50%. However, this 
decrease is compensated partially in k by a similar decrease in the standard deviations of eq. (P). 
As a result, the cut | ?? |< 1.75 decreases the absolute value of k by only 15%. 

If the average charged multiplicity (nt) is measured in the target region as a function of the 
multiplicity ric in the current region, the restriction | 7/ |< 1.75 decreases the (nt) by a factor of 2 
to 5, depending on ric. Therefore, for the present target-region acceptance, the average number of 
particles in the target region as a function of the multiplicity in the current region (or vice versa) 
cannot be reliably determined, which is why k is used, as described in subsection |2.3| . 



7 Uncorrected observables 



To understand how well the MC models describe simple observables related to the angular and 
current-target correlations, we compare the uncorrected data to the MC distributions after the 
detector simulation. Table |l| gives the numbers of events in each region of Q^, the average values 
of and the average charged multiplicity in the current region of the Breit frame for the data. 

Fig. ^ shows the probability, P„, of detecting n charged particles in the current region of the Breit 
frame for four regions. The difference between the data and Monte Carlo models is illustrated 
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by the x^/NDF in Table |^. Note that, although the angular correlations are sensitive to Pn-, they 
are also determined by the semi- inclusive two-particle densities pO[| . 

The average charged multiplicities in the current (target) hemisphere as a function of the multi- 
plicity in the opposite hemisphere are shown for > 10 GeV^ in Fig. ^. In both cases the average 
multiplicity in a hemisphere decreases as the number of particles in the opposite hemisphere in- 
creases, which is a signature of anti-correlations between current and target charged multiplicities. 

These comparisons show that ARIADNE gives the best description of the uncorrected observables. 
Therefore, this model is used to correct the correlations for detector effects. 



8 Correction procedure 

Due to the complexity of the measured variables, a bin-by-bin correction procedure is used. The 
correction factors C for each kinematic region in (Q^, x) are evaluated separately for each observable, 
A = f{e),r{e),K, 

/I gen 

C = ^^ (10) 

where A^'^'^ is calculated at the generator-level of ARIADNE and A*^^^ is that at the detector-level 
of this model. The corrected value for an observable is found by multiplying its measured value by 
the relevant correction factor. 

The correction factors are close to unity for all observables and vary smoothly for any given A. 
For f{e) and r(e), the correction factors are never larger than 1.3 and decrease with increasing 
e. For the current-target correlations, the correction factors never exceed 1.35. The sign of these 
correlations remains unchanged after the correction procedure. 



9 Statistical and systematic uncertainties 

The main sources of systematic uncertainties are: 



• Event reconstruction and selection. The systematic check was performed by varying the cuts 
on 7/ei yjB, S and the vertex position requirement (ye < 0.90, yjB > 0.05, 40 < 5 < 55 
GeV, -35 < 

■^vertex 45 cm). The contribution of these uncertainties for the angular and 
current-target correlations is typically 60% of the total systematic errors. 

• Track reconstruction and selection. The cuts were tightened: tracks should have transverse 
momenta larger than 200 MeV and | 77 |< 1.5. Tracks which reach at least the third CTD 
superlayer were used. These uncertainties are typically 30% of the total systematic errors. 

• The use of tracks not fitted to the primary vertex. 



The overall systematic uncertainty is determined by adding the uncertainties discussed above in 
quadrature. The error bars on the corrected data presented in Sect. 10 include both statistical and 
systematic errors, added in quadrature. The errors on the angular correlations are dominated by 
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the statistical errors (about 80 — 90% of the total error). For the current-target correlations, the 
major uncertainty comes from the systematic effects. 

No systematic uncertainty is attributed to the use of LEPTO or HERWIG in place of ARIADNE 
to determine the correction factors since these models do not adequately describe the raw data for 
several variables relevant to the correlation study. If the difference between the correction factors 
determined using ARIADNE and those using LEPTO or HERWIG is added, the systematic errors 
increase by about 50%. 



10 Results 



10.1 Two-Particle angular correlations 



The behaviour of the normalised particle density, ^"(e), (see eq. @) and correlation function, Y{e), 
(see eq. @), measured in the four regions, listed in Table |l], are shown in Figs. ^ and ^ for 
three values of 0. For these figures, the value of A = 0.15 GeV is chosen in order to be consistent 



with similar measurements made at LEP [12,13| 



The values of Y{e) increase with increasing e, reflecting a strong decrease of the number of particle 
pairs with increase of e. Such a trend is mostly determined by the single-particle distribution, rather 
than correlations between particles. In contrast, Y{£) is more sensitive to two-particle correlations. 
It rises with increasing e at large where the particle multiplicity is larger and the jet structure 
is more pronounced. This behaviour implies an increase in the strength of the correlations with 
decreasing angular separation between the two particles. At low Q^, however, this rise is not 
observed due to the small charged particle multiplicity in the current region. 

The Y{e) distributions of the MC models agree reasonably well with the data. The ARIADNE 
model describes the data in all regions. At low LEPTO slightly overestimates Y[e), while 
HERWIG underestimates it. The agreement is less satisfactory for Y{£). ARIADNE is quite 
successful in the description of low and medium Q^, but underestimates the data at high Q^. 
LEPTO fails to describe the data at < iqo GeV^. HERWIG gives a poor description of y(e) 
in all regions, except for the highest Q^. 

In the current version of ARIADNE 4.8, the suppression of available phase space for parton radiation 
due to the extended nature of the target remnant can also affect the current region of the Breit 



frame at high Q^. A recent (as yet unreleased) modification of ARIADNE |47|, which ensures 
that the whole of the available current-region phase space is used for gluon radiation, gives a good 
description of Y(e) at high (not shown). 

The results indicate that there are no changes in the Y{£) and Y{e) distributions at large £ as the 
angle Q is varied in the range from 45° to 90°. This observation supports the predicted scaling 
property of eqs. (Q) and (^). To compare the analytic calculations with the data, we choose O = 60°. 
Fig. P shows Y{£) and the asymptotic QCD predictions (Q) for fixed (dashed lines) and running 
(solid lines) strong coupling constant^ for three different values of A. For the calculations with a 
fixed coupling constant, the values A = 0.05,0.15,0.25 GeV for > 2000 GeV^ correspond to the 
effective as = 0.107,0.130,0.144 for the lowest order QCD relation between the as and A. The 

^ An evolution of the parton shower at a fixed energy scale is characterized by a running coupling constant as 
which reflects a change of energy scale as the parton shower evolves. 
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assumption of a fixed coupling constant is inconsistent with the data for all intervals, while 
the calculation with a running coupling constant shows reasonable agreement with the data at 
g2 > 2000 GeV^ and large A. 

A similar comparison of the analytic calculations for Y{e) with the data is shown in Fig. |^. As 
for Y(e\ discrepancies are observed at small . The agreement is better at high Q^; however, 
the calculations still overestimate the data at small e. This discrepancy is likely to be related to 
the neglect of energy- momentum conservation in the DLA, as recently discussed in 13|. This 



simplification makes the angular correlations in the analytic calculations more prominent for large 
angular separations, Q\2 (small e). 

Figs. |8| and ^ show our results for Y(e) and Y(e) at G = 45° and A = 0.15 GeV together with the 
DELPHI data at P = 45.5 GeV = 91 GeV) ll|. The ZEUS data are shown for > 2OOO 
GeV^, which corresponds approximately to P = {Q)/2 = 30.6 GeV. Despite the difference in 
energy, the results for Y{e) agree. Note that no significant energy dependence of Y{£) was found 



at LEP |T1. 



According to [T^, the behaviour of Y{£) at ^/s = 183 GeV is steeper than that at ^/s = 91 GeV. 
If there is a universality between the current region of DIS and a single hemisphere of e~'~e~, one 
might expect that e^p data at P = (Q) /2 = 30.6 GeV should exhibit a less strong rise than e+e" 
data at P = 45.5 GeV. However, the ZEUS data shown in Fig. ^ exhibit the opposite trend, which 
suggests differences between the angular correlations in the current region of DIS and a single 
hemisphere of e^e~ annihilation. The observed discrepancy might be related to different choices of 
the axis of the 0-cone: whereas this axis in e^e~ was determined by the sphericity axis, the virtual 
photon direction is used to determine the axis of the G-cone in DIS. In addition, the leading-order 
QCD effects discussed in Sect. |2| can lead to an uncertainty in our results. To investigate this 
issue, similar studies of the angular correlations were performed using only single-jet events. The 
single-jet pre-selection described in Sect. ^ leads to smaller values of Y{e) than the data shown in 
Fig. ^, but a small discrepancy with e^e~ still remains. 

Several further checks have been performed to understand the difference between the ZEUS and 
the DELPHI results for Y{£). Firstly, ln(PG/A) was used in the denominator of (|l|), as in the 
DELPHI analysis. Secondly, a small fraction of events at > 8000 GeV^ was rejected, so that 
no event has an initial parton at an energy higher than P = 45.5 GeV. The analysis was repeated 
by calculating e with P = {Q'^)/2, rather than using P = Q/2 for each individual event. In 
addition, HERWIG and the modified version of ARIADNE without the suppression effect in the 
current region were used. Both models describe Y{e) at > 2000 GeV^, and, therefore, they 
are better suited to correct the data in this region. For the checks discussed above, no significant 
changes in the ZEUS data were observed which can account for the discrepancy. 

Figs. ^ and|^ also show the analytic QCD results ||l8| at infinite energy and finite energy (P = 30.6 
GeV) separately for quark and gluon jets. The analytic prediction for gluon jets shows better 
agreement with the data than that for quark jets. However, it again fails to describe Y^e) at small 
e. A possible source of the discrepancy with the finite-energy quark prediction is the ratio P = 9/4 
of the mean parton multiplicity in gluon and quark jets used in the calculations for quark jets. If 
one goes beyond the DLA, the value of this ratio is smaller than 9/4. For example, in the next-to- 



leading log approximation, P ~ 1.8 |48|. A smaller value of R will bring the analytic prediction for 



quark jets closer to the gluon prediction [32|. 
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10.2 Current-Target correlations 



Figure shows the behaviour of the correlation coefficient k as a function of the average values of 
and x. The same bins in and x were used as in previous studies [||,|2|. Instead of calculating 
K in each bin, the bins are combined to increase the statistics. Table |2| gives the numbers of events 
and the uncorrected average values of and x. Fig. p!o| a shows the dependence of the correlations 



on Q , while Fig. ICb shows the x variation. Note that the average values of Q and x shown in 



the figure are corrected. 

Anti-correlations (k < 0) are observed for all values of x and Q^, as predicted by eq. @). The 
magnitude of k decreases with increasing (Q^) from 0.35 to 0.1. According to the analytic result 
of eq. @, the observed anti-correlations can be due to the 0{as) effects (QCDC and BGF). Their 
kinematics in the Breit frame can reduce the particle multiplicity in the current region and increase 



it in the target region. A Monte Carlo study |29| shows that the contribution from hadronisation 
is relatively small. 



The ZEUS result, shown in Fig. IC, gives a quantitative estimate of the correlations between the 
current and target region multiplicities in DIS at HERA, an effect which is quite different from 
that observed between the forward and backward hemisphers in e^e~ annihilation [^,p^25|. The 



kinematics of the 0{as) effects in the Breit frame which lead to these correlations are also a possible 
source of the disagreement between ZEUS and e^e~ data for the angular correlations. These effects 
can also lead to discrepancies between the data and the analytic QCD calculations which do not 
take into account the kinematics of the 0{as) processes in the Breit frame. It is noteworthy that 
the current region of DIS shows a different average multiplicity at small than a single hemisphere 
of e+e~ annihilation [@,^,@ due to the 0{as) processes. 

The magnitude of the anti-correlations increases with decreasing (x). According to the analytic 
result dll), this can be due to an increase of the fraction of events with one or two jets produced in 
the target region. This behaviour is driven by an increase of the gluon density inside the proton, 
leading to an increase of the Boson-Gluon Fusion rate. 

The ARIADNE model agrees well with the data. A reasonable description of dijet production 
in DIS by this model [49| is possibly responsible for this agreement. The LEPTO and HERWIG 



predictions show the same trend but do not reproduce the magnitude of the correlations which is 
likely to be related to a lower predicted dijet rate E^]. 



11 Conclusions 

The evolution of two-particle angular correlations with of DIS has been studied in the Breit 
frame. The data have been compared to the results of numerical simulations and to those of analytic 
QCD calculations in the DLA for partons, assuming the Local Parton-Hadron Duality hypothesis. 

The results on the angular correlations support the predicted scaling behaviour in the angle Q at 
large e. The scaling in energy is found to hold approximately for Y{e). For this variable, our results 
are in good agreement with those of DELPHI for e^e~, despite the difference in the energy. In 
contrast to Y{£), the two-particle correlation function Y{e) depends strongly on Q^: at low Q^, 
the two-particle correlations are suppressed, but they increase with e at high (Q^ > 100 GeV^). 
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The analytic calculations for Y{e) and Y{s) do not describe the data at low Q^. The results 
become closer to the data as increases. The asymptotic QCD predictions with a running 
coupling constant only describe the Y{s) at > 100 GeV^ for large values of the effective A. For 
the highest studied, the predictions reproduce the trend of Y{£) but overestimate the data at 
small e. 

The best description of the angular correlations is achieved by the ARIADNE model, although 
it fails to describe Y{e) at the highest studied. The prediction of ARIADNE shows smaller 
correlations in this region than the data due to the suppression of available current-region phase 
space by the remnant in the colour-dipolc model. LEPTO and HERWIG show small discrepancies 
for Y{e). For the Y{£) distribution, which is more sensitive to the two-particle correlations, both 
models fail to describe the data at low Q^. 

The current-target correlations in the Brcit frame arc found to be large and negative. They are thus 
very different from those measured in e~^e~ collisions, where small and positive forward-backward 
correlations have been observed. In DIS, the correlations increase with decreasing x. Such a 
behaviour is expected from an increase of the Boson-Gluon Fusion rate that leads to an increase 
of dijet production. ARIADNE agrees well with the current-target correlations. Neither LEPTO 
nor HERWIG reproduce the magnitude of the correlations. Insufficient dijet production in these 
models is a possible source of this failure. 

Many studies at HERA have found similarity between multiparticle production in the current 
region of e+p collisions and a single hemisphere of e'^e" annihilation. This paper demonstrates 
that correlations are a powerful tool for investigating this area. For high Q"^, the ZEUS results 
show differences between the two-particle angular correlations in DIS and e+e" annihilation; this 
is inconsistent with the universality of two-particle inclusive densities, which are expected to be 
sensitive to the kinematics of the first-order QCD processes producing final-state hadrons in both 
the current and target regions. 
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(GeV^) range 


(Q^) GeV^ 


(n) 


No. of events 


10-20 


14.4 ± 0.5 


0.96 ±0.04 


161421 


20-100 


40 ± 1 


1.5 ±0.1 


159453 


> 100 


320 ± 11 


2.9 ±0.3 


31461 


> 2000 


3737 ± 203 


3.9 ±0.6 


456 



Table 1: Numbers of selected events used to study the angular correlations in different regions, 
integrated over all x values. The corresponding uncorrected average values of and average 
charged multiplicity (n) in the current region of the Breit frame are also shown. The errors are the 
combined statistical and systematic uncertainties. 



(GeV^) range 


X range 


(Q") GeV^ 


(x) 


No. of events 


evolution with 


10-20 


0.6-2.4-10"^ 


14.2 ±0.6 


(1.22 ±0.01) - 10"^ 


79381 


20-40 


1.2 - 10 • 10^^ 


28.7 ±0.7 


(2.88 ±0.01) - 10"^ 


48287 


40-80 


1.2 - 10 - 10"^ 


55 ±2 


(3.82 ±0.01) - 10"^ 


34956 


80-160 


2.4- 10- 10~^ 


109 ±3 


(5.29 ±0.02) - 10"^ 


11331 


160-320 


2.4 - 50 - 10"^ 


218 ± 5 


(2.1 ±0.1) - 10"^ 


4350 


320-640 


10 - 50 - 10"^ 


416 ±9 


(2.4 ±0.2) - 10"^ 


1615 


640-1280 


10 - 50 - 10"^ 


899 ± 14 


(1.7 ±0.2) - 10"^ 


270 


evolution with x 


10-1280 


0.6 - 1.2 - 10"^ 


20.7 ±0.4 


(8.7 ±0.1) - lO"'' 


72039 


10-1280 


1.2-2.4-10"^ 


28 ±2 


(1.7 ±0.1) - 10"^ 


74621 


10-1280 


2.4-10-10"^ 


53 ±5 


(4.4 ±0.1) - 10-^ 


73825 



Table 2: Bins in and x used to study the current-target correlations. The average values of 
and x are shown without detector corrections. The errors are the combined statistical and 
systematic uncertainties. 



(GeV^) range 


XVNDF 


ARIADNE 4.8 


LEPTO 6.5 


HERWIG 5.9 


10-20 


25.4 


189.6 


73.8 


20-100 


27.7 


167.9 


71.1 


> 100 


4.5 


34.3 


32.9 


> 2000 


1.0 


1.3 


0.9 



Table 3: = (-^r?'^*'^ ~ -^^ti^^)^ /(^Data + ^Mc) P^'- degree of freedom for the uncorrected 

probability distribution P„ shown in Fig. 2|. 
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current 



target 




Figure 1: A schematic representation of the measurement of correlations in the Breit frame. 
The angular correlations in the current region are measured between any two charged particles 
separated by 9i2- Long-range correlations in the full phase space are measured between current 
(nc) and target (rtt) multiplicities. 
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Figure 2: Probability distribution P„ for detecting n charged particles in the current region of 
the Breit frame for different ranges of Q^. The uncorrected data are compared to MC predictions 
after the detector simulation. The errors on the data are the statistical uncertainties. The shaded 
bands show the statistical uncertainties on the MC predictions. Table ^ gives the x^/NDF of the 
MC predictions. 
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Figure 3: Average charged multiplicities in the current (target) hemisphere as a function of the 
multiplicity in the opposite hemisphere for > 10 GeV^. The uncorrected data are compared to 
MC predictions after the detector simulation. The statistical uncertainties on the data are typically 
smaller than the size of the symbols. The shaded bands show the statistical uncertainties on the 
MC predictions. 
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Figure 4: The normalised particle density Y{e) in the current region for four ranges and three 
different values compared to MC predictions. The rescaling is performed using A = 0.15 GeV. 
The inner bars on the corrected data show the statistical uncertainties. The full error bars include 
the systematic uncertainties, which are typically negligible compared to the statistical errors. 
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Figure 5: The correlation function Y{e) in the current region for four ranges and three different 
G values compared to MC predictions. The rcscaling is performed using A = 0.15 GeV. The inner 
error bars on the corrected data show the statistical uncertainties and the full error bars include 
the systematic uncertainties. 
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Figure 6: The normalised particle density Y{e) in the current region for = 60° compared to 
the analytic QCD predictions (^) for different effective A values. Solid (dashed) lines show the 
predictions at asymptotic energy for a running (fixed) coupling constant. 
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Figure 7: The correlation function Y{e) in the current region for = 60° compared to the analytic 
QCD predictions (|5|) for different effective A values. Solid (dashed) lines show the predictions at 
asymptotic energy for a running (fixed) coupling constant. 
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Figure 8: The normalised particle density Y{e) for > 2OOO GeV^ (P = (Q)/2 = 30.6 GeV) 
compared to the analytic QCD predictions and DELPHI results for P = 45.5 GeV. For both 
experiments = 45° and A = 0.15 GeV. The solid line shows the DLA predictions at asymptotic 
energy for a running coupling constant. The dashed (dotted) line corresponds to the prediction at 
finite energy (P = 30.6 GeV) for quark (gluon) jets. 
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Figure 9: The correlation function Y{e) for > 2000 GeV^ (P = {Q)/2 = 30.6 GeV) compared 
to the analytic QCD predictions and DELPHI results for P = 45.5 GeV. For both experiments 
= 45° and A = 0.15 GeV. The solid line shows the DLA predictions at asymptotic energy for a 
running coupling constant. The dashed (dotted) line corresponds to the prediction at finite energy 
(P = 30.6 GeV) for quark (gluon) jets. 
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Figure 10: (a) represents the evolution of the coefficient of correlations k with predominant 
variation in for corrected data and MC predictions; (b) shows the same quantity where pre- 
dominantly X varies. The corrected values of (<3^) and (x) are indicated for each plot. The inner 
error bars on the data show the statistical uncertainties. The full error bars include the systematic 
uncertainties. 
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